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Abstract 

We use a functional approach to evaluate the Casimir free energy for 
a self-interacting scalar field in d + 1 dimensions, satisfying Dirichlet 
boundary conditions on two parallel planes. 

When the interaction is turned off, exact results for the free energy in 
some particular cases may be found, as well as low and high temper- 
ature expansions based on a duality relation that involves the inverse 
temperature /3 and the distance between the mirrors, a. 
For the interacting theory, we derive and implement two different ap- 
proaches. The first one is a perturbative expansion built with a ther- 
mal propagator that satisfies Dirichlet boundary conditions on the 
mirrors. The second approach uses the exact finite-temperature gen- 
erating functional as a starting point. In this sense, it allows one to 
include, for example, non-perturbative thermal corrections into the 
Casimir calculation, in a controlled way. 

We present results for calculations performed using those two ap- 
proaches. 

1 Introduction 



Casimir and related effects, where quantum effects depend upon the existence 
of boundary conditions for a quantum field, have been extensively studied [I]. 
Many different points of view and approaches to this kind of problem have 
been followed, under quite different sets of assumptions regarding the system; 
i.e., its intrinsic properties, and the conditions under which one wants to 
evaluate the Casimir force. 

In particular, there has been much interest in studying the Casimir effect at 
a finite temperature (T > 0), a study that can be undertaken at different 



1 



levels, distinguished by the way of taking into account all the possible T > 
effects. One could, for example, include thermal effects in the description of 
the matter on the mirrors (a very active research topic [2j), or rather for the 
physical description of the vacuum field [3], or for both. 

In this article, we shall assume perfect (at any temperature) mirrors, with 
thermal effects restricted to the vacuum field. 

It has been noted that thermal and Casimir-like effects share some similar- 
ities; this should hardly be surprising, since both may be regarded, essentially, 
as finite-size effects, the former in the Euclidean time interval [0, [3] (with pe- 
riodicity/antiperiodicity conditions), and the latter for a spacial coordinate 
with Dirichlet or Neumann conditions. 

Even though the nature of the boundary conditions is different, it should 
be expected that, when both effects are present, interesting relations ('dual- 
ities') will arise between the dependence on the inverse temperaure (3 = 
and the distance between the mirrors. 

In this article, we present an approach to the calculation of the Casimir 
free energy for a scalar field at finite temperature in o? + 1 dimensions, sub- 
ject to Dirichlet boundary conditions, including self-interactions. This ap- 
proach is based on the use of a path-integral formulation, whereby the d+1- 
dimensional problem is mapped to a dimensionally reduced one, with fields 
living on the boundaries. 

For all the cases considered we analyze the high and low temperature ex- 
pansions, computing also perturbative corrections to the Casimir free energy 
due to the non-interacting fields. 

This work is organized as follows: in section [2] we present our approach, 
introducing conventions and definitions. In section [3] we deal with the calcu- 
lations and the corresponding results for a free real scalar field. Interactions 
are introduced in [Hand, in [5], we present our conclusions. 



2 The method 

The main object we shall be interested in is the free energy F{P, a) for a real 
scalar field ip in d + 1 spacetime dimensions, which is subject to Dirichlet 
boundary conditions on two plane mirrors, separated by a distance a along 
the direction corresponding to the Xd coordinate. 

In terms of the corresponding partition function Z{P, a), F{l3, a) is given 

by: 

F{P,a) = -I In Z{P, a), (1) 
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while for Z{[3, a) we shall use its standard functional integral representation: 

Z(/?,a) = j [V^]e-'^^^ , (2) 

where S is the Euclidean action at T > 0; i.e., with an imaginary time 
variable restricted to the [0, /?] interval [4] and periodic boundary conditions 
for the fields with respect to this coordinate are implicitly assumed. 

In this work, we consider an action S such that S = Sq + Sj, with the 
free part of the action, 5*0, is given by: 

Soiv) = I jy^J d''x{d^ifd^^ + m\'') , (3) 
while the interaction part, S"/, is given by: 

Si{^) = ^jy^J d'x[v{T,^)]\ (4) 

We have used brackets in ['Py^j to indicate that the path-integral measure in- 
cludes only those field configurations satisfying Dirichlet boundary conditions 
on the locii of the mirrors, which correspond to two parallel planes. Accord- 
ingly, we use a coordinate system such that, if the Euclidean coordinates are 
denoted hy x = (xq, Xi, . . . , Xd) {xq = r), the mirrors then correspond to the 
regions: Xd = and Xd = a, and will be parametrized as follows: 

X = (a;||,0) , X = {x\\,a) , (5) 

where x\\ = (r,X||) = (r,xi, . . .,Xd-i). 
Then we have: 

[Dip] = Vif 6[ip{T,:s^i\,0)] 5[(/?(r,X|i,a)] , (6) 

where the (5's with braketed arguments are understood in the functional sense; 
for example: 

5[(/?(r,X||,a)] = JJ 5(^(r,X||,a)) , (7) 

-r,x|| 

while the ones on the right hand side are ordinary ones. Besides, we as- 
sume that the (/^-dependent) factor that comes from the integration over the 
canonical momentum has been absorbed into the definition of Vip, so that 
performing the integral over ip does indeed reproduce the partition function 
(without any missing factors) . 
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Following [5], the functional (5-functions are exponentiated by means of 
two auxiliary fields, ^ind ^2{x\\), living in d spacetime dimensions: 

Z(P,a) = y"l?^I?eiPe2 e-^(^) + ^^^'^'^-^''(^)^(^), (8) 

where we introduced the singular current: 

Jp{x) = + S{xd - a)^2{xi\) , (9) 

whose support is the region occupied by the mirrors. 

Since we shall treat interactions in a perturbative approach, it is conve- 
nient to deal first with the free theory, as a necessary starting point to include 
the interactions afterwards. 



3 Free theory 



When the theory is free (A = 0), 5* = 5*0, and the integral over ip becomes 
Gaussian. The resulting partition function, denoted by Z^^\i3,a) is then: 



(10) 



where Z^^\(3) is the (free) partition function in the absence of the mirrors, 
and: 

1 



SAO 



where a, 6 = 1, 2, and 



d x\\d y\\^a{x\\)Vtah{x\\,y\\)ii{y\\) 



A(a;||,0;y||,0) A(a;||, 0; a) 
A(x||,a;y||,0) A(a;||, a; yy, a) 



(11) 



(12) 



where A is the free, imaginary-time, propagator. It may be written explicitly 
as follows: 

A(r., x; Ty,y) = -^Y.j (^^^^'^^ -^)+^'^(''-^) A(.;„, k) , (13) 
with A(a;„,k||) = {ul + £\^)y\ 

We have used the notation £^(k||) = y''k| + m^, and Un = ^ denotes the 

Matsubara frequencies. A is the inverse, in the space of r-periodic functions, 
of = i-d'^ + m^). 
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Expression (fTO|) allows one to extract from the free energy the term that 
would correspond to a free field in the absence of mirrors, F^^\ plus another 
contribution, which we shall denote by 

F<»'((3,a) = F'°>(/3) + F^'HAa). (14) 

Fr(fta) = -iln[^^]. (15) 

which clearly contains the Casimir effect information (including thermal cor- 
rections). However, it also carries information that is usually unwanted, 
associated to the ever-present self-energy of the mirrors. Indeed, we see that 
it includes a divergent contribution to the free energy, which can be neatly 
identified, for example, by noting that, when a oo: 

(16) 

where is the contribution corresponding to one mirror: 



1)^2 e"^-^''''^"^^^^!!^^^^!!''''^"'''^^''^^"^ . (17) 



We then identify 

f!^\P) = -^InKH/^)] (18) 

as the free energy term that measures a mirror's self-interaction. It is, of 
course, independent of a. Extracting also this contribution, we have the 
following decomposition for F^^^ 

F('\p,a) = fW(/?) + 2F^\P) + F^'\p,a) (19) 

where Fc^\(3, a) has a vanishing limit when a ^ oo. 

We can produce a more explicit expression for the interesting term Fc^\ 
starting from its defining properties above. Indeed, we have: 

f!:'\M = -1 In I A , (20) 

^ ^ P U(0)(/3)[Z(?)(/3)]^J ' ^ ^ 

which, by integrating out the auxiliary fields, may be written in terms of the 
matrix kernel Q: 

F!:'\(3,a) = ^Trlnl] - ^Trlnfio, (21) 



where Q^c = ^U^oo and the trace is over both spacetime coordinates and 
a, b indices. 

Assuming that an UV regularization is introduced in order to make sense 
of each one of the traces above, and using 'reg' to denote UV regularized 
objects, we see that: 



1 
1 

2P 



Trlnfi 



reg 



Trln {Q^'Q) 



1 

reg 



Tr In Qr 



J reg 



(22) 



As we shall see, the trace on the second line is finite, and has a finite limit 
when the regulator is removed. Using some algebra we may reduce the trace 
to one where only continuous indices appear: 



2/3 



Tr In Qr 



where a 'reduced' kernel flc which is given by: 

^c(a;||,i/||) = -yii) - T(x|j 

with: 



(23) 



(24) 



d'^z\id'^w\\d\i[Qii] ^(x|j, ^11)^12(2:11, 



(25) 



has been introduced. 

Finally note that, due to translation invariance along the mirrors, the free 
energy shall be proportial to Vp, the area of the mirrors. Then, to absorbe 
this divergence we shall rather consider the corresponding free energy density, 
obtained by dividing the extensive quantity by 1^. In particular, 



JP-i°)(/5,a) = lim 



{(3, a) 



(26) 



where F^^^^ on the right hand side is evaluated for a system with a finite 
parallel volume. Moreover, taking advantage of the translation invariance 
along the parallel coordinates, we may use a Fourier transformation to write: 



+00 



2(3 J (27r)^- 



1 y: \nnrih)^ 



(27) 
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where 

^c^\\h\) = 1 -T(")(k||), (28) 

with: ^ _ _ _ _ 

TW(k|j) ^ [fiS^(k||)]-^fii?(k||)[fi(?(k||)]-i(](^(k||). (29) 

where the tilde denotes Fourier transformation in both r and xy. Note the 
appearance of the reciprocals of the matrix elements of d^''^^ (not to be con- 
fused with the matrix elements of the inverse of that matrix). The explicit 
form of the objects entering in (l29l ) may be obtained from ( fT2l) : 



nS)(k||) = i]g)(k|| 



1 



-a^a;2+f2(k||) 

fil?(k,|)=f](^Hk||) = -^===f=^, (30) 
ni^Hkii) = 1 - e-'V-'+^'C'ii) . (31) 



so that: 



It is then evident that the sum and the integral in (|27[) converge, since 
the integrand falls off exponentially for large values of the momenta/indices. 
This behaviour should be expected, since we have constructed this object 
subtracting explicitly the would-be self-energy parts, the possible source of 
UV divergences. 

Thus, the main result of the previous calculations is a (finite) expression 
for the free energy, which can be written as follows: 

") = i / 7^ E to [l - e-^V^^ . (32) 

^ ' n=— CO 

It is worth noting that the expression above has been obtained subtracting 
the free energy corresponding to a situation where the mirrors are separated 
by an infinite distance; on the other hand, we note that, when the /3 — oo 
limit is taken, the sum is replaced by an integral, and we obtain: 

)(a) . (oo, a) = ^ / 1^ In (l - e'^'^^) , (33) 

which is the Casimir energy per unit area. 

With this in mind, one may introduce yet another quantity; the tempera- 
ture dependent part of the free energy, whose area density we will denote by 
Tt^\(3,a)^ and, from the discussion above, it is given by: 

(/5, a) = (/3, a) - J-f (oo, a) . (34) 
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By using a rather simple rescaling it the integral in the second term, we may 
write: 

^ ' n=— OO 

/oo 
du In [1 - e-fV(2-)^+(/'^(''ii))^] I , (35) 
■oo 

which yields the Casimir free energy as a difference between a series and an 
integral, although the for the thermal part, since the sum over energy modes 
has already been (implicitly) performed. 

Before evaluating the free energy for different numbers of spacetime di- 
mensions, we explore below some consequences of a duality between (3 and 
a. 



3.1 Duality 

The dual role played by /? and a in the free energy may be understood, for 
example, by attempting to calculate that object by following an alternative 
approach, based on the knowledge of the exact energies of the field modes, 
emerging from the existence of Dirichlet boundary conditions. 
The energies uji of the stationary modes are: 



u;,(k||) = y— ^ + kj + m2 , leN. (36) 

Since each one of these modes behaves as a harmonic oscillator degree of 
freedom, its free energy / [tf;;(k||)] has the following form: 

^ oo 

/h(k||)] = -ilnfX;^-'^'"'^'"^^'^^^^ 

^ N=0 

= l^^,(k||) + iln[l - e-^^'^'^^O] . (37) 

Now we consider Tf'\j3, a), the temperature dependent part of the free energy 
density, obtained by summing the second term in the expression above over 
all the degrees of freedom, and dividing by the parallel area: 
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This may also be written as follows: 



e 



We now recall that we are considering the free energy for one and the same 
system, albeit with different normalizations: subtraction at a ^ 00 in (l32l) . 
and at /3 — s> 00 in ( l39l) . Then we may write: 

^W(/3,a) - lim (/3,a) = ^r(/5,a) - lim , (40) 

or: 

^f(/3,a) - ^W(a) = J-f (/3,a) - ^f)(/3) , (41) 

where is an a-independent object, which corresponds to the free 

energy of a free field. 

Thus, if we are going to be concerned with temperature and a-dependent 
quantities, for example, to study the temperature dependence of the Casimir 
force, we only need derivatives with respect to /3 and a of the expression 
above, and we see that: 

K^HAa)] =7JL[-^i°HAa)], (42) 



dadp^ " 'J dp da 
or: 

(/?, a) = (/?, a) ^ ^(0) ^) ^ (43) 

where the tildes denote subtraction of any term which has a vanishing mixed 
second partial derivative. 

On the other hand, from the explicit form of the free energy density as 
sum over modes (Matsubara or Casimir), we find the identity: 

[/3^f](2a,/3/2) = (/5, a) , (44) 

which, combined with fH3l) yields a duality between a and /?: 

^W(2a,/3/2) = A^W(/5,a). (45) 

To proceed, we make the simplifying assumption that m = 0, and transform 
variables in the integral over ky to have a dimensionless integral, obtaining 

P^^'\p,a) = (^f-^g{j,d), (46) 
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where g{'-f, d) is given by 

^?(7, d) =CdY. r 1^ - e-^^/^) , (47) 

n 

where 7 = Ana/ [3 and is a constant factor that depends solely on the 
dimension d. 

Now, using the duality formula fllSl) . we obtain the an interesting relation 
involving g: 

Y-'9il,d) = a'-'gia^d), (48) 

where a = Pit /a (or a = (27r)^/7). Note that this relation has immediate 
relevance to relate the low and high temperature regimes. Indeed, writing 
the expression above more explicitly, we see that: 

,«) 

This duality relation had been pointed out, for the d = 3 case, by Balian 
and Duplantier [6]. In the remainder of this section, we apply the previous 
results and study some particular properties of the free energy corresponding 
the the free scalar field, firstly for the d = 1 case, which we single out since 
it can be exactly solved, and then for c? > 1. 

3.2 d=l 

The free energy (neglecting an irrelevant constant) is given by the expression: 
J-W(/3,a) = i-J]ln(l-e-2'^l-l), (50) 



which may be written as an infinite product 

+00 



•^f = ^ In 11(1- g^"), (51) 

n=l 

where q = e~^™/^, and | g | <1 for T > 0. Using standard properties of 
ellipthic functions, we get a more explicit results for the free energy: 

•^i°nA«) = ^ + ^ln[r/(2a//5z)], (52) 

where ri{z) is Dedekind's eta function. Although it has been obtained for 
T > 0, one can verify that it also yields the proper results for T — > 0, namely. 

Besides, the ?7-function satisfies the property: r]{l/z) = \fizr]{z) which, 
in this context, is tantamount to the duality relation. 
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3.3 d>l 



Now, we exploit the duality relation to extract the high and low temperature 
behavior of the Casimir energy, for d> 1. In a high temperature expansion, 
formula (i47l ) can be expanded for small e'"' . The infinite-temperature limit 
corresponds to 7 ^ oo, thus n = yields the leading contribution in this 
expansion, with the n = 1, 2, . . . terms producing higher-order corrections. 
The explicit form of the leading high temperature term is: 

P^)Wa) L_imm p^o (53) 

whereas aplying the duality formula ( !45l) a — > /3/2, we have the first order 
contribution in the low temperature regime, 

Tf\2a,PI2) ~ (p ^ (54) 

which is a independent. The sub-leading contributions {n = 1,2, ...) include 
more involved integrals; nevertheless, when d = 3 and 7 — > 00, it reduces to 
terms like e~"^. Thus, the most significant contribution for T — 00 is: 

Therefore, the duality formula implies that, when T — ^ 0, 

•^f (A«)~-^e--^/^ (56) 
Collecting these results for d = 3, we see that, for high temperatures: T — > 00 

while, for low temperatures the corresponding behaviour is: 

360 720 



^f(2a,/5/2) . 



1440a3 



l + ^^C(3) + ^^e-"'3/'^ 



(58) 



4 Interacting theory 

Let us now consider the theory that results from turning on the quartic self- 
interaction term, studying the corresponding corrections to the free energy. 

This can be done in (at least) two different ways. We shall consider first 
the would-be more straightforward approach, which amounts to expanding 
the interaction term in powers of the coupling constant first, postponing the 
integration over the scalar field until the end of the calculation. We call it 
'perturbative' since it is closer in spirit to standard perturbation theory. 
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4.1 Perturbative approach 

We have to evaluate: 

= Z(°)(/5,a) (e-^^(^)), (59) 

where we have used (. . .) to denote functional averaging with the free (Mat- 
subara) action and the scalar-field integration measure satisfying Dirichlet 
boundary conditions. Namely, 

Thus, the different terms in the perturbative expansion are obtained by 
expanding in powers of the coupling constant, A. For example, the first order 
term, J^^^\ may be written as follows: 

J-«(/?,a) = A| d'+'xMx)]'), (61) 

which may be expressed (via Wick's theorem), in terms of the average of just 
two fields. Indeed, the generating functional for these averages, (/3, a), is 
simply: 

Zi^{P,a) = j [p<^]e-^o(^)+/'^'*+'-'^(-M-) . (62) 

We again introduce the auxiliary fields to impose the periodicity constraints 
on the measure, so that: 



Integrating out the scalar field ^p, and using a simplified notation for the 
integrals, we obtain: 

Zo{P,a) = Zo{p,a) j exp -]- j ^„(x||)fi«/3(x||, y||)^/3(i/||) 

+ if Um)La{x\\) + I [ J{x)Aix,y)J{y)\ , (64) 

^X|| ^ J x,y 

where: La{x\\) = A(x|| , Oq,; y) J('y), with = (a — l)a, a = 1,2. Integrat- 
ing now the auxiliary fields, we see that: 

Z^iP, a) = ZoiP, a) ^^>y J(-)Gi.,y)Jiy) ^ ^g^^ 
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where G, the (free) thermal correlation function in the presence of the mir- 
rors, may be written more explicitly as follows: 



G{x;y) = {Lp{x)ip{y)) = A{x;y) - / M(x||, x^; xj,, a^) 

X [^"^]a/3(x||,y||)A(y||,a^;y||,yrf)| , (66) 

and this is the building block for all the perturbative corrections. 

It is quite straightforward to check that the propagator above does ver- 
ify the Dirichlet boundary conditions, when each one of its arguments ap- 
proaches a mirror. For example: 

lim G{x;y) = A(x||,a^;i/) - / \n.fa{xi\]Xu)[n~%p{x\\,yu)A{y[,,a^-,yi\,yd) 



= A{x\\,a^]y) - / 6^i36{x\\,y\^)A{y\^,a(^]y\\,yd) 

= A{x\i,aY,y) - A{xi\,aY,y) = 0, (67) 

and analogously for the second argument. 

The thermal correlation function G{x; y) above is composed of two contri- 
butions: one is the usual free thermal propagator in the absence of boundary 
conditions and the other come from reflexions on the surface boundaries. In 
this way we write 

Gix;y) = A{x-y)-M{x-y), (68) 
where M(x; y) is given by: 

M{x;y) = ^Y. j --^)+^'^ii(''ii--^ii)M(^„,k|,;x,,7/,) , (69) 



with M(w„, k||; Xrf, i/d) equal to: 

Q-{\^d\+\yd\)E _ ^-{\xi-a\+\ya\+a)E _ ^-{\xj\ + \ya-a\+a)E _|_ ^-{\xa-a\ + \yd-a\)E 

2E{1 - e-2'^^) 



(70) 



where E = + ky + m^. 

We apply now the previous approach to the calculation of the first-order 
correction to the free energy (in fact, to the force) and the self-energy func- 
tion: 
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4.1.1 First-order correction to the free energy 

We want to evaluate thermal corrections to Casimir energy in the interacting 
theory, in particular, to understand how its divergences should be dealt with. 
To automatically subtract a-independent contributions, we work with the 
derivative of the free energy with respect to a. Obviously, this is a force, 
and it has the same amount of information as an energy which has been 
subtracted to avoid a-independent infinities. Indeed, it can be integrated 
over any finite range of distances to obtain de energy difi'erence. 
Using Fcas to denote these derivatives, we have: 

F = + F^^) (71) 

^ cas cas cas ' \ ' 

Where the term Fell is the first-order correction to free Casimir force Fc^. 
There is a term where only A(a:, y) appears; this can be renormalized as usual 
in finite-temperature field theory, by a zero-temperature subtraction plus 
the inclusion of the first-order (temperature-dependent) mass counterterm 
contribution [4]. 

Then, keeping temperature and a dependent terms only, we see that 

_ I (70) 

-^cas ~ -^cas, AM ^ -^cas,MM ' V^y' 

where: 

^ii',.M = ~A-(.,.)/.-.^^^ (73) 

^i„M =|;/<i-.M(x.^)^^t^. (74) 

Using the notation /a,m and /a/ a/ for the corresponding area densities (we 
omit super and subscripts): 

= -^A (0) ^ J J ^ dx, , (7o) 



and 



and 



n,l 

X M(k||,w;„,Xd) ^ , (76) 

and the first-order correction to the Casimir force will be given by: 

/a,m + /m,m • (77) 



^(1) 



A. 



d-l 
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The term M(ty„, ky; Xf^) can be obtained from (I70|) : after changing vari- 
ables from X — > a; + I (to obtain a symmetrized form), we perform an 
integration over d + 1-dimensional spacetime, obtaining: 



/a,m 



The T-dependent function A-'"(0) is given by 

A^(0) = 



- csch^ (ai?fc 



coth(a-E'fc) 
2E, 



(78) 



(79) 



where nB{P,uj) is the Bose distribution function, with u = i/p^ + m^. 

/a, A/ may be expressed in yet another way, using the definition of the 
Bose distribution function but with P replaced by 2a while u is replaced by 

= \/^n '^'^l^ "^sil^-i^) is replaced by 



nk{2a,Ek) = 



1 



Then 



/a. 



M 



Tp 



A^(o)$: 



(2vr; 



d-l 



Eh 



2ank (1 + Uk] 



(80) 



(81) 



which is clearly a convergent quantity. 

For the remaining term, we proceed in a similar way, where now the p 
represent {uji, p\\) 



(82) 



where 



[rir. 



Tlr 



1)1 ^ 
1 



UkU. 



v 



El Ek/2a 

rik Ek rip - Ep Uk 



" ' 2 J \EkEp Ek{El - Ej) J 2Ek{Ek + Ep 
which is also convergent. 



(83) 
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4.1.2 First-order correction to the self-energy 

We conclude the application of this method with the calculation of the tad- 
pole diagram contribution to the self-energy. There is here a surface diver- 
gence due to the M{x; x) contribution: 

n = ^(A(x,x) - M(x,x)) . (84) 

As before, it is convenient to separate from the self-energy the contribution 
which is present even when the mirrors are infinitely distant; namely, 

n = Hfree + ^mir , (85) 

where Hfree = jA{x,x) and Umir = —jM{x,x) which contain the contribu- 
tion to the selfenergy coming from the Dirichlet boundary conditions. The 
term Hjree has UV and IR divergences, which are usually analyzed in stan- 
dard finite temperature calculations. 

On the other hand, the Umir term has no UV divergences, although it 
is IR divergent in low dimensions. Moreover, it presents surface divergences 
when the Xd coordinate tends to = or Xd = a. We will classify these 
divergences according to the dimension of the theory. 

When d = 1 and m = 0, we obtain, 

^mir = -^Y. 2E(1 - e-2«^) ' ^^^^ 

where E = \ujn\- Then we see that 

if s < 
if < s < 1 

e-7(^-l)n if S > 1 



An ^-^ n 

n=l 



cogfe((s-l/2)7n)-e-T"/^ 
sinh(jn/2) 



where 7 = Aira/P and s = x/a. The sum for s < and s > 1 can be 
performed. For s < 0, the selfenergy is proportional to ln(l — exp(7s)), which 
diverges for s 0, regardless of the temperature. 

For (i > 2, it is convenient to separate the self-energy part which into two 
pieces, one coming from n = 0, and the other coming from the remaining 
modes: 

TT — TT"=0 _|_ TT^T^O (orr\ 
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and 



( g2sp 



rmrr 



if s < 



-—— / dp - X <( 



coshii2s-l)p)-e-P ifO<S<l 
smn(p) 



'2{s-l)p 



if s > 1 



n: 



X 



2-Kd 5^ 

n=l 



if s < 



X < 



Y^p2„i sinh("inp/2) 



if s > 1 



t i^o(n7(s - 1)) 

Again, the zero mode contribution nj^~?^ has an IR divergent behaviour 
for low momenta, while there is a divergence on the mirrors. For instance, 
for s < the selfenergy is proportional to log{—2s). 

The n^^°^ term can be analized in the high and low temperature limits, 
i.e. 7 S> 1 and 7 -C 1, respectively. The divergence of the function Kq{x) is 
logarithmic in a; = 0, thus surface divergences also are logarithmic. Moreover, 
we see that > n^t°. 

For d > 2, the IR divergences are absent, and the surface diverces appear 
explicity. The two contributions are given by: 

^ if s < 



n: 



n=0 



Ar(f-i) 

p2'^+^n^a'^'^ 



n: 



C{d - 2, 1 - s) + C(rf - 2, s) - 2C{d - 2) if < s < 1 

if s > 1 

if s < 

cos/i((s-i)7np)-e ^"^ 



/ ^ynsp 



d-3 
2 < 



sm/i{7n£) 
g-7n(s-l)p 



if < s < 1 
if s > 1 



The integral for s < or s > 1 can be performed 

' (_i)f-if^izitZ!!l i 

,00 \ -^J 4_i 



n: 



A 



/3t (2a)f 



X < 



n=l 



d -1 



if s < 



if s > 1 
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In the case of ci = 3 dimension we have a simple expression 



1 



Ma(3 I 1 



(s_l)(e7(s-i)„i) 



if s < 



if s > 1 



By the above results we conclude that the surface divergences of the n^~°^ 
term are polynomial: ~ s"'^'^"^-', whereas the term I^^^r which has the sum 
of non-zero modes presents a more severe divergence proportional to s^*^'^"^-'. 
This has just been shown explicity for the c? = 3 case. 

In figures 1, 2, 3 and 4, we plot the self-energy for 1, 2, and 3 dimensions, 
for different values of the parameters. 



4.2 Non-perturbative approach 

Let us conclude this section by considering a second, alternative procedure 
to the one just explained. It amounts to integrating out, albeit formally, the 
scalar field before the perturbative expansion. Indeed, introducing, from the 
very beginning, the auxiliary fields used to impose the Dirichlet conditions, 
we see that the partition function becomes 



where W denotes the generating functional of connected correlation func- 
tions, at finite temperature, and with an unconstrained (no Dirichlet condi- 
tions) scalar field integration measure: 



VLpe^^^'^^^J •^^'^''^^'^> . (89) 

Of course, the arbitrary current J in the definition above must be replaced by 
Jp, which does depend on the auxiliary fields and on a, the distance between 
the two mirrors. Besides, we assume that, in the course of evaluating W, 
a renormalization procedure has been used to make sense of the possible 
infinites, in the usual way. 

To proceed, we recall that W[J] does have a functional expansion: 

W[J] = Wo + ^j d^+^x j d''+'yW2{x,y)J{x)J{y) + ... (90) 

where Wk correponds to the connected fc-point correlation functions. Odd 
terms are, for the quartic perturbation, absent from the expansion. Then one 
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may invoque some approximation that allows one to truncate the functional 
expansion. Of course, it cannot be a naive perturbative expansion in the 
coupling constant; one should rather use, for example, a mean field or large- 
expansion. Then the leading term will be just the quadratic one: 

Z{f3,a) ~ Z{f3) X Z,{(3,a) (91) 

where Z{(3) = is the thermal partition function in the absence of mirrors, 
while 

Z,{(3,a) = j V^,V^2 e-'^S'i'^'^S'i'^'ym[^^y)M^)My) ^ (92) 

depends on W2{x,y), the full renormalized thermal propagator. Using the 
explicit form of Jp, we see that the integral is a Gaussian: 

Z,(/5,a) = j Pei^^6e-^««) (93) 

where 



S,{0 = I I d%d\U^)Ql,{x^\,y^\My\\) , (94) 



(95) 



and 

Qqt^ ■y)= \ D{xii,0;yii,0) D(x||, 0; yy, a) 

where D is the full imaginary-time propagator. It may be written quite 
generally as follows: 

Dir., x; w) = Wj T^e^^"^^^ 5(^„, k) , (96) 

n 

with _D(co'„,k||) a generally complicated function of its arguments. However, 
we note that some non perturbative corrections do produce a simple result. 
For example, considering the IR resummed version of the massless scalar 
field, yields an expression of the form: 

D(a;„,k) = [ul + kl + U^{T)]-' (97) 

where n^(T) is the thermal mass. For example, the first two non-trivial 
contributions correspond to a term which is linear in A plus a non-analytic 
term: 

n,(T) = ^[l-3(^)i + ^^0^ (98) 



19 



Upon insertion of this expression, we see that the corresponding contribution 
to the Casimir free energy becomes: 



1 



J^c{P,a) = — 



-2a /a;2+k2+^^(T)■ 



2/5 J {271^- 

where the a — > oo contribution has already been subtracted. 



(99) 



5 Conclusions 

We have obtained exact results for the free energy in a Casimir system in 1 + 1 
dimensions and well as low and high temperature expansions for li > 1, based 
on a duality relation between inverse temperatures and distances between 
mirrors, in d + 1 spacetime dimensions. 

For the interacting theory, we have derived two different approaches. In 
the first approach, perturbative expansion has been used to obtain the first 
order corrections to thermal Casimir force, showing that it is finite once the 
standar renormalization of thermal field theory in the absence of mirrors is 
performed. 

On the other hand, we also studied the selfenergies, showing that they 
have surface divergences in their coordinate dependence, and that those di- 
vergences are of a polynomial type, with a degree which depends on the 
dimension. 

Finally, we have shown that a variation in the order used to integrate the 
auxiliary fields yields a different method, whereby the one particle irreducible 
functions naturally appear. 
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Figure Captions 



Figure 1: Self-energy in 1 + 1 dimensions, for three values of 7 = 
Ana/ (3: upper curve: 7 = 10; middle curve: 7 = 1 and lower curve: 
7 = 0.1. 

Figure 2: Self-energy in 2 + 1 dimensions, for three values of 7 = 
Ana/ (3: upper curve: 7 = 10; middle curve: 7 = 1 and lower curve: 
7 = 0.1. 

Figure 3: Self-energy in 3 + 1 dimensions, zero mode contribution. 

Figure 4: Self-energy in 3 + 1 dimensions, for three values of 7 = 
Ana/ (3: upper curve: 7 = 10; middle curve: 7 = 1 and lower curve: 
7 = 0.1. 
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